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Abstract. Simplified model for current-voltage characteristics of weak links is 
suggested. It is based on an approach considering the multiple Andreev reflection 
in metallic Josephson junction. The model allows to calculate current-voltage 
characteristics of the superconductor - normal metal - superconductor junctions with 
different thicknesses of normal layer at different temperatures. A hysteretic peculiarity 
of V(I) dependence is described as result of the negative differential resistance. The 
current-voltage characteristic of high-Tc composite YBCO -l-BaPbOs were computed. 
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1. Introduction 

Superconductor - normal metal - superconductor (SNS) junctions have the current- 
voltage characteristics (CVCs) with rich peculiarities. Given certain parameters of 
junction, CVCs of SNS junctions demonstrate the excess current, the subharmonic gap 
structure and the negative differential resistance at low bias voltage. The region of 
negative differential resistance corresponds to the hysteresis of voltage in the bias current 
measurements. The SNS junctions with nonlinear CVCs are promising for different 
applications, e.g. low- noise mixers in submillimetre-wave region switcher 
nanologic circuits 

Description of CVCs of SNS junctions was subject of many articles and there 
were recognized the key role of multiple Andreev reflections |S1 13 El 13 CHI- The 
main features of CVCs enumerated above are successfully described by Kiimmel - 
Gunsenheimer - Nicolsky theory (KGN) j7j. KGN theory is applicable for the thick 
and clean weak links, where the normal metal layer N has the thickness 2a larger than 
the coherence length of superconductor, the inelastic mean free path / larger than 2a. 
A simplified model in frame of KGN theory was developed by L.A.A. Pereira and R. 
Nicolsky [11 . This simple model is relevant for the weak links with thin superconducting 
banks S. The contribution of scattering states is omitted in the model |llj. 

Pereira - Nicolsky and KGN model were applied earlier to describe experimental 
CVCs of various weak links [121 ESI UH UH UHl El- Experience of applications 
demonstrates that oversimplified Pereira - Nicolsky model gives only qualitative 
description. We suggest a new simple modification of KGN theory. It is shown that 
the CVCs of SNS junctions can be computed without all the complex Ansatz of KGN 
theory. We hope this will lead to more extensive using of KGN based approach to the 
calculation of weak link characteristics. 

2. Current- voltage characteristics 

2.1. Model 

Let us consider a voltage-biased SNS junction with a constant electric field which is in 
negative z direction perpendicular to the NS interfaces and exists in the N layer only 
(Figure P). The normal layer has the thickness 2a. The thickness of superconducting 
bank is D — a >> 2a. 

Dynamics of quasiparticles in the SNS junction was considered in the work [Zj, 
where the time dependent Bogoliubov - de Gennes equations are solved for the wave 
packets of nonequilibrium electrons and holes. The expression for dissipative current 
density in SNS junction was resulted in [7]. Calculations of density of states [THj and 
probability of finding of the quasiparticles [7] in the N region are requested preliminary 
to calculate the current. The expression for current through SNS junction with thick 
superconducting banks {D — a » 2a) can be written as following [7j: 
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where gr{E) is the two dimensional density of states, Pn{E) is the probabihty of finding 
of the quasiparticles with the energy E in the N region, m* is the effective mass of 
electron, / is the inelastic mean free path and Rn is the resistance of the N region with 
the thickness 2a, A is the value of energy gap of superconductor at the temperature 
T, kzF is the z component of Fermi wave vector of quasiparticles, n is the number 
of Andreev reflections which quasiparticles undergo before they move out of the pair 
potential well. 

The probability Pn{E) of finding of quasiparticles with energy E in the N region 
is given by Eq.(2.19) of [3: 

with the penetration depth A = ^ ^a'^-e ^ ior E < A, \ < D — a and \ = D — a 
otherwise. For the quasiparticles from the scattering states Pn{E) = 2a/2D. Let us 
accept for the sake of simplicity A >> a. Therefore Pn{E) = 2a/2\ for the bound 
states. 

2.2. Density of states 

The density of states is found from 
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where A is the normal layer area, kzF,r defines the value of kzF for which Er = E. 

The energy spectrum Erik-^p) consists of the spatially quantized bound states and 
the quasicontinuum scattering states. The energy eigenvalue equation for the spatially 
quantized bound Andreev states [7j is transcendental; it calculated numerically only 
Figure El 
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where r = 0,1,2,. . . 

Let us simplify Eq.Q. The expansion of arccos(£'r/A) in to Taylor series 
(7r/2 — Er/A +. . . ) up to the second term and the subsequent expressing of Erik^p) 
are executed. Then we inserted the correcting multiplier C for best fitting of Eq.(P)). 
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Figure 2. Energy of the bound Andreev state with r — 2a ~ 5000 A; A = 0.57 
meV; kp = 1.62 A'^. 1) Eq.©, C = 0; 2) Eq.®, C = 1; 3) the exact solution of 
Eq.Q; 4) Eq.©, C = 7r/2(l - am* A/h^kp); 5) Eq.©, C = tt/2. 

If C = then the spectrum of Pereira - Nicolsky model is reproduced (curve 1, 
Figure El). R. Kiimmel used Eq.(ini) with C = tt/2 (THl for approximated calculation of 
the energy spectrum (curve 5, Figure El). 

The density of bound states follows from ^ and (jSJ 

^ (e, 
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For quasiparticles from the quasicontinuum states, the energy spectrum is 
approximated by the continuous BCS spectrum of homogeneous superconductor [TllTHj: 
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For SNS junction with thick superconducting banks, the effective energy gap A* 
equals A. Then the density of quasicontinuum scattering states is 

A2m* ^ E 

g{E) = — — -fc^D ^ (8) 
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Figure 3. Density of states g{E) of SNS junction witli tliick superconducting banks 
resulted by JSI (dotted line), g{E) calculated by Eq.Q with C = 7r/2 (dash line), 
g{E) calculated by Eq.® with C = 7r/2(l - am* A/fi^kF) (solid line). D = 70000 A; 
2a = 5000 A; = 3.77 K; kp = 1.62 A-\ 

It is reasonable to choose the variable multiplier C from the best fitting of the energy 
spectrum and the density of states. We suggest C = vr/2(l — am* A/h^kp) for C > 1 
and C = 1 otherwise. Such choice of C provides a good agreement of Eq.© (curve 
4, Figure El) with the numerical solution of Eq.(j3|) for different relation of a,m* , A, kp. 
The coincidence of density of states resulted with g{E) by [THJ is satisfactory (Figure Ej). 



2.3. Current density 

The current density of quasiparticles from bound states is resulted with ((H) and (jHl): 
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After substitution of (jSj) to ((H) and exclusion of small term, we get the current 
density of quasiparticles from quasicontinuum states: 
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where Ei = —A + neV for —A + neV > A and Ei = A otherwise. 

The current densities and (fTUI) include the vohage dependence only within the 
integral limits. 

If eV » ksT, A the integral in (fTUj) can be transformed and the excess current 
density is resulted: 

^»W-^^.=Ata„h(^)exp(-^) (11) 

This excess current density is the same as one obtained in KGN theory (Eq.(4.12) 
in HI). 

Note that jbsiV) dependence does not change practically if the second summation 
in © is interrupted at r = 0. Therefore we can write the expression for total current 
density as follows 
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C = 7r/2(l - anfA/h^kp) for C > 1 and C = 1 otherwise; £"1 = -A + neV for 
— A + neV > A and £"1 = A otherwise. 

This simplified model allows to calculate the CVCs of weak links with thick 
D ~ a » 2a superconducting banks. It operates for different thicknesses of normal 
layer 2a < I and different temperatures lower Tc. 



3. Comparison with experimental current-voltage characteristics 

The model is successfully applied for the description of CVCs of SNS junctions pUj . 
The steep rise of current density at low voltage, the arches of subharmonic gap structure 
(SGS), and the excess current are reproduced (Figure in the computed CVCs (fT^ . 
The multiple Andreev reflection is the main reason for these peculiarities (Zj. Position 
of the (n-l)th arch of SGS is between Vn and Vn+i, where Vn = 2A/{n — l)e, so the 
largest 1th arch is between A and 2A. Small peaks on the arches near Vn are caused by 
subgap peak on g{E). 

The number of allowed Andreev reflections decreases with increasing of bias voltage 
jZj. The current density due to Andreev reflections decreases correspondingly. In that 
time the ohmic current density increases as well as voltage. Region of the negative 
differential resistance appears in the CVC if the decreasing of Andreev current density is 
stronger than the increasing of the ohmic current density. The measurement of current- 
biased CVC demonstrates the hysteretic voltage jumps (Figure Ej) instead the region of 
negative differential resistance. 
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Figure 4. Voltage-biased current-voltage characteristic of SNS junction with steep 
rise of current (a), region of negative differential resistance (b), excess current density 
jex (c). Tc = 3.77 K, Ao 0.57 meV, hp = 1.62 A-\ 2a = 5000 A, / = 15a, T = CIT^. 
Arrows (d) and (e) display the hysteretic jumps of voltage in current-biased CVC 

To prove the model we firstly [20] described the CVCs of tin microbridges PHI^ I^ 
in different temperatures. The hysteretic peculiarity is absent in the CVCs. Satisfactory 
agreement of the calculated curves and the experimental data is achieved pUj for the 
known parameters of Sn. 

The application of model is possible to describe the CVCs of combination of weak 
links e.g. networks and contacts connected in series. The networks of weak links, which 
realize in the polycrystalline high-Tc superconductors, have the composite CVCs. These 
CVCs are a superposition of individual CVCs of the single weak links that constitutes 
the network. Fitting parameters for current and voltage should be used in the model 
to account a straining of CVC along / and V axes |^ . The weak links connected 
in series (SNSNS...) can be realized in the break junctions and wires with phase slip 
centers. In this case, the straining of CVCs along V axis is accounted by the formula 
for series of weak links with dispersion of parameters I15| ITH]. 

We used the developed model to compute the CVC of composite 92.5% YBCO -|- 
7.5% BaPbOa ^Sj that is the network of weak links. We get the satisfactory agreement 
(Figure EI) for the parameters of investigated composite (Tc = 93.5 K, Aq = 17.5 meV, 
kp = 0.65 2a = 50 A, / = 9a, m* = Am^, T = 4.2 K). The scale of right axis in 
units eV/A demonstrates the straining of CVC due to the superposition of individual 
CVCs of single weak links. This description of experimental CVC is more correct than 
earlier one by Pereira-Nicolsky model. 

4. Conclusion 

The simplified model for calculation of current-voltage characteristics of SNS junctions 
was developed. The KGN approach j7j was changed to be more convenient for 
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Figure 5. Current- voltage characteristic of composite YBCO + BaPbOa. Experiment 
[Oj (points) and calculations (solid lines). 

description of experimental CVCs of weak links with thick superconducting banks. 
The model operates for different thicknesses of normal layer 2a < / and different 
temperatures lower T^. The frequent observed peculiarities (steep rise of current, arches 
of subharmonic gap structure, negative differential resistance, excess current) in CVCs 
of SNS junctions are interpreted to be produced by the multiple Andreev reflections. 
The hysteretic peculiarity is described as result of the negative differential resistance. 

The model was applied to compute the hysteretic current- voltage characteristics of 
high-Tc composite YBCO -fBaPbOg. 
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